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The relativistic six-quark equations are constructed in the framework of the dispersion relation 
technique. The relativistic six-quark amplitudes of dibaryons including the light u, d and heavy c, 
b quarks are calculated. The approximate solutions of these equations using the method based on 
the extraction of leading singularities of the heavy hexaquark amplitudes are obtained. The poles of 
these amplitudes determine the masses of charmed and bottom dibaryons with the isospins / = 0, 
1, 2 and the spin- parities = 0+, 1+, 2+. 
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. I. INTRODUCTION. 

■ 

' In 1977, JafFe T| studied the color-magnetic interaction of the one-gluon-exchange potential in the multiquark 
^\ \ system and found that the most attractive channel is the flavor singlet with quark content u^(Ps^. The same symmetry 
analysis of the chiral boson exchange potential leads to the similar result [2]. 

The iJ-particle, A^fi-state and di-ft may be strong interaction stable. Up to now, these three interesting candidates 
^ ^1 of dibaryons are still not found or confirmed by experiments. It seems that one should go beyond these candidates 

1 ' and should search for the possible candidates in a wider region, expecially the systems with heavy quarks, in terms 

of a more reliable model. 

^ There were a number of theoretical predictions by using various models: the quark cluster model 0, the quark- 

delocation model ^, the chiral SU{3) quark model 0, the fiavor SU{3) skyrmion model Lomon predicted a 
deuteronlike dibaryon resonance using R-matrix theory [9|. By employing the chiral SU{3) quark model Zhang and 
Yu studied ilfl and Sil states [l^, El- Lee and Yasui discuss the stable multiquark states containing charm and 
: bottom quark [l^ . 

' In a series of papers [T3l - [l7l | a method has been developed which is convenient for analyzing relativistic three-hadron 
, systems. The physics of the three-hadron system can be described by means of a pair interaction between the particles, 
ps) ■ There are three isobar channels, each of which consists of a two-particle isobar and the third particle. The presence 
Q,^ of the isobar representation together with the condition of unitarity in the pair energies and of analyticity leads to a 

■ system of integral equations in a single variable. Their solution makes it possible to describe the interaction of the 
^-H ' produced particles in three-hadron systems. 

y—i I In our papers [l8l - l20j relativistic generalization of the three-body Faddeev equations was obtained in the form of 
^ ■ dispersion relations in the pair energy of two interacting quarks. The mass spectrum of S'-wave baryons including w, d, 
• ^ ,: s quarks was calculated by a method based on isolating the leading singularities in the amplitude. We searched for the 

■ approximate solution of integral three-quark equations by taking into account two-particle and triangle singularities, 
^ , and defined all the smooth functions of the subenergy variables (as compared with the singular part of the amplitude) 
- - I in the middle point of the physical region of Dalitz-plot, then the problem was reduced to the one of solving a system 

of simple algebraic equations. 

In the present paper the relativistic six-quark equations are found in the framework of coupled-channel formalism. 
We use only planar diagrams; the other diagrams due to the rules of 1/Nc expansion 21-23] are neglected. The 
six-quark amplitudes of dibaryons with two heavy quarks are calculated. The poles of these amplitudes determine 
the masses of heavy dibaryons. 

In Sec. II, the six-quark amplitudes of hexaquarks are constructed. The dynamical mixing between the subampli- 
tudes of dibaryons are considered. The relativistic six-quark equations are constructed in the form of the dispersion 
relation over the two-body subenergy. The approximate solutions of these equations using the method based on the 
extraction of leading singularities of the amplitude are obtained. Sec. Ill is devoted to the calculation results for the 
dibaryon mass spectra (Tables HI HIl) . In conclusion, the status of the considered model is discussed. 
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II. SIX-QUARK AMPLITUDES OF THE HEXAQUARKS WITH THE TWO HEAVY QUARKS. 

We derive the relativistic six-quark equations in the framework of the dispersion relation technique. We use 
only planar diagrams; the other diagrams due to the rules of l/iV^ expansion [2l| - |23| are neglected. The current 
generates a six-quark system. The correct equations for the amplitude are obtained by taking into account all 
possible subamplitudes. Then one should represent a six-particle amplitude as a sum of 15 subamplitudes: 



(1) 



i<3 
i.3 = l 



This defines the division of the diagrams into groups according to the certain pair interaction of particles. The 
total amplitude can be represented graphically as a sum of diagrams. We need to consider only one group of diagrams 
and the amplitude corresponding to them, for example A12. We shall consider the derivation of the relativistic 
generalization of the Faddeev-Yakubovsky approach. In our case, the low-lying dibaryons with the two heavy quarks 
are considered. We take into account the pairwise interaction of all six quarks in the hexaquark. 

For instance, we consider the state EcSc with the isospin 1 = 2 and the spin-parity ~ 0+ {uuc uuc). The 
set of diagrams associated with the amplitude A12 can further be broken down into seven groups corresponding to 
subamplitudes: Af", A}", Af\ ^r""o"<^^ ^o^'o^" ^ ^i^.^.^icc^ 

The system of graphical equations (see for example equation for the amplitude " for the state Sc^c with the 
isospin 1 = 2 and the spin-parity = 0+ {uuc uuc)) is determined by the subamplitudes using the self-consistent 
method. The coefficients are determined by the permutation of quarks. 

In order to represent the subamplitudes A\ , A\ , A^ ^ A\ ^ , A\ ° , A^ ^ , A\ ^ ^ in the form of a 
dispersion relation, it is necessary to define the amplitude of qq, qQ and QQ interactions. We use the results of our 



relativistic quark model 2J| and write down the pair quark amplitudes in the form: 



an{s^k) = 5-7 T , (2) 

1 - Bn{Stk) 



A 

^2 /„/ 



J TT S,. - Sik 



(mi+mt,)2 



^ y/ jsik - (™^ + mk)'^){sik - (m^ - mfc)^) 

•^ik 

The coefficients a{n,J^), /3{n,J^) and d{n,J^) are given in Table Hill Here n = 1 coresponds to qq, qQ and 
QQ-pairs with = 0+, n = 2 corresponds to qq, qQ and QQ with = 1+. 



The coupled integral equations correspond to Fig. 1 can be described similar to [25|. 

Then we can go from the integration of the cosine of the angles dzi to the integration over the subenergies. 

Let us extract two- and three-particle singularities in the amplitudes Al , Al , Af , A\ ^ , A\ ° , A§ , 
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N _ <-i=2 ■'12345, •'1234, ^12, ^34;-uiu„ ^ai27^1""V^347 

(.,S12345,S1234,S12,.34) - [I - B,..{s,Ml - B,..[s,,)] ' 

.l^-'O"'^^ \ ar"°"'(s, S12345, S1234, Si2, S34)Bi>...(si2)Bo"c(s34) 

^2 Gs,Sl2345,Sl234,Si2,S34) = t- rry- , (9) 

.0"'=0"'^/ \ 02"'""' (s, 512345, Sl234,Si2,S34)So"<:(si2)So"c(s34) 

A, {S, .12345, .1234, .12, .34) = [1 - So„. (.,2)] [1 ^ So-^ (.34)] ' ^'^^ 

s ar"^""^"(s,Sl2345,Sl2,S34,S56)Sii.^(si2)Si>.,.(s34)Bicc(s56) 

^3 (s, .12345, .12, .34,556) = 



[1 - Si.. (.12)] [1 - Si.. (.34)] [1 - Blcc(s56)] 



(11) 



We used the classification of singularities, which was proposed in paper \2d\. Using this classification, one defines 
the reduced amplitudes ai, 02, as as well as the _B-functions in the middle point of physical region of Dalitz-plot at 
the point sq. 

Such choice of point .0 allows us to replace integral equations (EcSc, 1 = 2, = 0+) by the algebraic equations 

(m - &■ 

af" = A + 4q}""/i(1""1"") +4af"/i(r"0"=) +2a^""^""/2(r"l""l"") +8a^""°"'/2(l""l""0"=) 

+ 2af''°'"/2(l""0"'=0"'=) (12) 

= A + 8ar^i(l"0"") + 12ar°"'/2(l""0""0"") (13) 

af = A + 3a}""/i(0"'=l"") +ar^i(0"'l'') +4af'/i(0"'=0"'=) + 6ai""°"''/2(0"n™0"'=) 

+ 3af' °'"'/2(0"^0"'=0"=) (14) 

^i-i"" ^ A + 4Q}""/3(l""l""r'") +8af 74(1""1""0"=) + 16af """77(l""r'"r'"0"=) 

+ af '°"'(4/5(l™l""0"^0"'=) +8/6(l""l""0"'=0"'=)) (15) 

= A + a}"" (2 /3(1""0""1"") +/4(0""1""1"")) +af' (2 /3(1""0""0"") + 2 /4(r"0""0"')) 

+ o:^""""*^ (2 /g(i""o"'=l""0'"=) + 2 /g(i""o"'0"n"") + 2 /^(i™o'"=l™0"'=) + 2 /7(i™o"^0"'=l"") 

+ 2/7(0"^l""l""0"'=)) +af'°"'(/5(0"=l""0"'=0"^) + 2/6(l""0"=0"^0"'=) + 2/7(0''^l""0"'=0"=)) 

(16) 

= A + a}""(/3(0"^0"n"") +4/4(0"^0"'=r")) +ar^3(0'"'0"'=l'='=) +aP(2/3(0"'=0"^0"=) 

+ 4/4(0'''=0''^0"'=)) + 2ai""i"76(0"'=0"n''"r") + """'(4/5(0"'=0''n"''0"'=) 

+ 4/6(0'"=0'"=l™0"'=) + 4/7(0"'=0"=r'"0"'=) + 4/7(0"'=0"'=0'"=l™)) + a^"'°"'(2/6(0"'=0"'=0'"=0"'=) 

+ 4/7(0'"=0'"=0'"'=0'"=)) + 2af "i""i"/8(0"'=0'"=l™r=r") (17) 

^i""i""i" ^ ^ ^ 4 Q,}""/g(i""i"«i<='=i™) + 8 af 79(i™'r'=i™o"") + 16 ap°"7io(i""r"i"'r"o"') 
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+ 8af'""'/io(l""l'='=r"0''^0'"=), (18) 

where are the current constants. We used the functions Ii, I2, I3, I4, h, le, I7, Is, I9, Iio similar to the paper [27[. 
The solutions of the system of equations are considered as: 



where zeros of D(s) determinants define the masses of bound states of dibaryons. 



III. CALCULATION RESULTS. 



The model in question take into account the hexaquarks with the two heavy quarks qqqqQQ, q — u,d, Q — c,b: 
uuuucc, uuudcc, uuddcc, uuuubb, uuudbb, uuddbb. 

The quark masses of the model are niq = 495 MeV, rric — 1655 MeV and mf, = 4840 MeV. 

The experimental data are absent, therefore we use the dimensionless parameters, which are similar to the previous 
paper [23 ■ It allows us to calculate the mass spectra of qqqqQQ states. We use the gluon coupling constants go = 0.653 
(diquark = 0"*") and gi — 0.292 (diquark — 1+), cutoff parameter A = 11. We consider the Aqc,cc — 8.52, 
which are determined by M = 5250 MeF (/J = 22 E^S*, the threshold is 5290 MeF. In the case of &-quarks 

the AqbM = 7.35 is determined by M = 11620 MeF (/J = 22 SfcS*, SJ^^), the threshold is 11660 AfeF. 

We have calculated the heavy dibaryon masses with isospin / = 0, 1, 2 and spin-parity — 0+, 1+, 2+, which 
are given in the Tables U and |IT1 The relativistic six-body approach possesses the dynamical mixing and allows us 
to calculate the contributions of the subamplitudes to the hexaquark amplitude (Tables |IV] - IVip . The calculated 
dibaryon subamplitudes A2 present the main contributions to the hexaquark amplitude (about 70 percents). We use 
only two new parameters for the calculation of 23 qqqqcc states and 19 qqqqbb states. 

The lowest mass for the qqqqcc states is M = 4364 MeV , for the qqqqbb states is M = 8670 MeV . 

In quark models, which describe rather well the masses and static properties of hadrons, the masses of the quarks 
usually have the similar values for the spectra of light and heavy hadrons. However, this is achieved at the expence 
of some difference in the characteristic of the confinement potential. It should be borne in mind that for a fixed 
hadron mass the masses of the constituent quarks which enter into the composition of the hadron will become smaller 
when the slope of the confinement potential increases or its radius decreases. Therefore, conversely, we can change 
the masses of the constituent quarks when going from the spectrum of light to the heavy hadrons, while keeping the 
characteristic of the confinement potential unchanged. We can effectively take into account the contribution of the 
confinement potential in obtaining the spectrum of heavy hadrons. We neglect with the mass distinction of u and d 
quarks. The estimation of the theoretical error on the heavy dibaryons masses is 1 MeV . This result was obtained 
by the choice of model parameters. 



IV. CONCLUSIONS. 



In a strongly bound systems, which include the light quarks, where p/m ~ 1, the approximation of nonrelativistic 
kinematics and dynamics is not justified. In our paper, the relativistic description of six-particles amplitudes of heavy 
dibaryons with the two heavy quarks is considered. We take into account the u, d, c, b quarks. Our model is confined 
to the quark- antiquark pair production on account of the phase space restriction. Here vriq and toq the "mass" of the 
constituent quark. Therefore the production of new quark-antiquark pair is absent for the low-lying hadrons. 

Hadronic molecules are loosely states of hadrons, whose inter-hadron distances are larger than the quark confinement 
size. 

The heavy analogue of H dibaryon (AcAc) does not exist though its potential is attractive 28;]. Oka et al. believe 
that the future studies may specify the binding energy of such a molecule state. The binding energy is sensitive to 
the cutoff parameter [2^ . Our calculation allows us to obtain the Ac Ac molecule bound state (Table H]) . There exist 
a loosely bound state with a small binding energy Eb = 54 MeV . In the case of ScSc dibaryon we obtain the bound 
state, but the binding energy is equal to Eb = 540 MeV . 

The similar results are obtained for the other heavy dibaryons with the two heavy quarks (Tables HI - HI]) . For 
the AfeAb system does not exist a loosely bound state with a small binding energy. The heavy dibaryons and heavy 
baryon-antibaryons may be produced at LHC. 
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TABLE I: S-wave charmed dibaryon masses. Parameters of model: cutoff A = 11.0 and Aqc.cc = 8.52, gluon coupling constants 
go = 0.653 and gi = 0.292. Quark masses rUq = 495 MeV and mc = 1655 MeV. 
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TABLE II: S-wave bottom dibaryon masses. Parameters of model: cutoff A = 11.0 and Aqb,bb = 7.35, gluon coupling constants 
go = 0.653 and gi = 0.292. Quark masses = 495 MeV and rub = 4840 Mel/. 
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TABLE 111: Vertex functions and Chew-Mandelstam coefficients. 
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TABLE IV: IJ = 00 SqEq, E^S^. Parameters of model: cutoff A = 11.0, A^ccc = 8.52, Aqi^u = 7.35, gluon coupling 
constants go = 0.653 and gi = 0.292. Quark masses niq = 495 MeV, rric = 1655 MeV and mb = 4840 MeV. 
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TABLE V: 7 J = 00 EqAq, AqAq. Parameters of model: cutoff A = 11.0, Aqc.cc = 8.52, Agb^bi, = 7.35, gluon coupling constants 
go = 0.653 and gi = 0.292. Quark masses rUq = 495 MeV, rUc = 1655 MeV and mb = 4840 Me V^. 
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TABLE VL IJ = 00 NEqq, ASqq. Parameters of model: cutoff A = 11.0, Aqc.cc = 8.52, Aq6,66 = 7.35, gluon coupling 
constants go = 0.653 and gi = 0.292. Quark masses niq = 495 MeV, rUc = 1655 MeV and mb = 4840 MeV. 
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